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Abstract 

Background: The consumption of a vertebrate blood meal by adult female mosquitoes is necessary for their reproduction, 
but it also presents significant physiological challenges to mosquito osmoregulation and metabolism. The renal 
(Malpighian) tubules of mosquitoes play critical roles in the initial processing of the blood meal by excreting excess water 
and salts that are ingested. However, it is unclear how the tubules contribute to the metabolism and excretion of wastes 
(e.g., heme, ammonia) produced during the digestion of blood. 

Methodology/Principal Findings: Here we used RNA-Seq to examine global changes in transcript expression in the 
Malpighian tubules of the highly-invasive Asian tiger mosquito Aedes albopictus during the first 24 h after consuming a 
blood meal. We found progressive, global changes in the transcriptome of the Malpighian tubules isolated from mosquitoes 
at 3 h, 12 h, and 24 h after a blood meal. Notably, a DAVID functional cluster analysis of the differentially-expressed 
transcripts revealed 1) a down-regulation of transcripts associated with oxidative metabolism, active transport, and mRNA 
translation, and 2) an up-regulation of transcripts associated with antioxidants and detoxification, proteolytic activity, 
amino-acid metabolism, and cytoskeletal dynamics. 

Conclusions/Significance: The results suggest that blood feeding elicits a functional transition of the epithelium from one 
specializing in active transepithelial fluid secretion (e.g., diuresis) to one specializing in detoxification and metabolic waste 
excretion. Our findings provide the first insights into the putative roles of mosquito Malpighian tubules in the chronic 
processing of blood meals. 
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Introduction 

The Asian tiger mosquito Aedes albopictus is considered one of tlie 
most invasive mosquito species in the world; since 1979 it has 
spread to over 28 countries outside of its native range in Asia and 
Southeast Asia, aided by the international trade of used 
automobile tires [1,2]. Within the United States, the mosquito 
has spread to at least 36 states and models of its potential for range 
expansion in the northeastern United States within the next few 
decades are alarming [3]. Moreover, this species is a known or 
suspected vector of several medically important arboviruses, 
including chikungimya, dengxie, eastern equine encephalitis, La 
Crosse, West Nile, and yellow fever [4]. Thus, A. albopictus is an 
emerging threat to global health for which effective control 
measures need to be developed. 

Historically, mosquitoes have been controlled through the use of 
insecticides that target the nervous system (e.g., carbamates, 
organophosphates, organochlorines, and pyrethroids). However, 
resistance to these control agents is limiting their efficacy. In 



particular, the yellow fever mosquito Aedes aegypti exhibits high 
levels of resistance to insecticides in certain parts of the world, and 
there is concern that A. albopictus will soon develop such resistance 
[5] . Thus, it is important to identify new control agents that target 
novel physiological systems in mosquitoes to help combat the 
emerging threat of insecticide resistance. 

A recent study by our group demonstrated that the renal 
excretory system (Malpighian tubules) of mosquitoes represents a 
valuable new physiological target for insecticides [6]. The 
Malpighian tubules produce urine via transepithelial fluid 
secretion, which is mediated by the coordinated actions of a V- 
type H -ATPase along with several ion transporters, ion channels, 
and water channels [7] . In adult female A. aegypti mosquitoes, the 
Malpighian tubules play an especially important role in the post- 
prandial diuresis when the mosquito excretes urine during and 
after the engorgement of vertebrate blood [8]. The diuresis lasts 
for up to two hours after feeding and excretes a significant fraction 
of the ingested Na"*", K""", CI , and water from the blood [8] . Once 
this diuresis ends, the female mosquito will continue to digest and 
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Author Summary 

The Asian tiger mosquito Aedes albopictus is a vector of 
several nnedically-important arboviruses and one of the 
most invasive mosquito species in the world. Existing 
control measures for mosquitoes are presently being 
challenged by the emergence of resistance to insecticides 
that target the nervous system. Thus, it is necessary to 
identify novel physiological targets to guide the develop- 
ment of new Insecticides. We recently demonstrated that 
the 'kidneys' (Malpighian tubules) of mosquitoes offer a 
valuable, new physiological target for insecticides. How- 
ever, our understanding of how this tissue contributes to 
the chronic metabolic processing of blood meals by 
mosquitoes is limited. Here we characterize the changes 
in transcript expression that occur in the Malpighian 
tubules of adult female A. albopictus with the goal of 
identifying key molecular pathways that may reveal 
valuable targets for insecticide development. We find 
dramatic changes in transcript accumulation in Malpighian 
tubules, which 1) provide new insights into the potential 
functional roles of Malpighian tubules after a blood meal, 
and 2) reveal new potential molecular pathways and 
targets to guide the development of new insecticides that 
would disrupt the renal functions of mosquitoes. 



metabolize the blood meal over the next two days to nourish the 
development of her eggs. The physiological importance of the 
Malpighian tubules during this time is unknown, but they 
presumably play a critical role in excreting excess nitrogenous 
wastes and other metabolites that are generated during the 
processing of the protein-rich meal [9] , especially within the first 
24 hours when ~75-90% of the ingested protein is digested 
[10,11]. 

Other groups have documented the effects of ingesting blood on 
the transcriptomes of adult female A. aegypti and Anopheles gamhiae 
mosquitoes [12,13,14,15,16], including more focused studies on 
how blood feeding influences tissue-specific transcriptomes in the 
antennae, fat body, midgut, and salivary glands of these species 
[17,18,19,20]. However, no previous studies have examined the 
effects of blood-feeding on the transcriptome of mosquito 
Malpighian tubules. 

The goal of the present study was to characterize the global 
changes in transcript expression that occur in the Malpighian 
tubules otA. albopictus during the first 24 h after female mosquitoes 
consume a blood meal (using RNA-Seq), with the aim of 
identifying key metabolic pathways and transcripts that are 
activated or suppressed in the renal tubules during the processing 
of the blood meal. We found that blood feeding elicits dramatic, 
time-dependent changes to the Malpighian-tubule transcriptome 
of A. albopictus. A functional cluster analysis of the differentially- 
expressed transcripts revealed a potential functional transition of 
the tubule epithelium after blood feeding from one specializing in 
active transepithelial fluid secretion to one specializing in 
detoxification and metabolic waste excretion. 

Methods 

Mosquito colony 

A. albopictus eggs were obtained from the Malaria Research and 
Reference Reagent Resource Center (MR4) as part of BEI 
Resources Repositoiy, NIAID, NIH (ALBOPICTUS, MRA-804, 
deposited by Sandra Allan). Eggs were raised to adults using a 
protocol similar to that described for A. aegypti [21] with the 



exception that larvae were fed pulverized TetraMin flakes (Melle, 
Germany). Adult females between 5 to 10 days post-eclosion were 
used for the present study. 

Design of the study 

The experimental design consisted of two treatments, blood fed 
(BF) and non-blood fed (NBF) females at three different time 
points. In brief, the BF mosquitoes were fed on heparinized rabbit 
blood for 30 minutes (see details below) and collected at 3 h, 12 h, 
or 24 h after feeding. These time points occur after the post- 
prandial diuresis, which ends within 2 h after feeding [8]. 
Moreover, one or more of these time points has been commonly 
used in other studies examining the effects of blood feeding on 
gene expression in mosquitoes [13,14,15,16,17,18,19,20,22]. NBF 
females were only offered a 10% sucrose solution and dissected at 
similar time points to serve as controls. Each treatment/ time point 
was replicated three times using females from different cohorts 
(i.e., 3 biological rephcates per time point). 

Blood feeding and isolation of Malpighian tubules 

For each experimental treatment, 90 adult female mosquitoes 
were transferred to two small 32 oz. cages (45 females per cage) 
without access to a sucrose solution for 24 h prior to offering them 
blood or sucrose. To one cage, a membrane feeder (Hemotek, 
Blackburn, UK) was used to feed the mosquitoes warmed blood 
(37°C), which consisted of heparinized rabbit blood (purchased 
from Hemostat Laboratories, Dixon, CA) supplemented with 
adenosine 5 '-triphosphoric acid (disodium salt; Sigma, St. Louis, 
MO) at a concentration of 0.01 g/ml. A solution of 10% lactic 
acid was applied to the membrane surface as an attractant. 
Females were given access to the blood for a period of 30 min 
before the feeder was removed from the cage. In the other cage, 
the mosquitoes were given access to cotton balls soaked with 10% 
sucrose for 30 min. At 3 h, 12 h, or 24 h after removing the feeder 
or cotton balls, the cage was refrigerated on ice to immobilize the 
mosquitoes. 

Before dissecting the mosquitoes that were offered blood, their 
abdomens were visually examined to confirm their engorgement. 
The alimentary canal of each mosquito was then extracted by 
tugging on the last segment of the abdomen with fine forceps 
under Ringer solution. The Ringer solution consisted of (in mM): 
150 NaCl, 3.4 KCl, 1.7 CaClz, 1.8 NaHCOg, 1.0 MgClg, 5 
glucose, and 25 HEPES (pH 7.1). The Malpighian tubules were 
isolated from their attachment to the alimentary canal and 
immediately immersed in 50 |.lL of TRIzol Reagent (Life 
Technologies, Carlsbad, BA) in a sterile 1.5 ml microcentrifuge 
tube on ice. A total of 200 Malpighian tubules (from 40 females) 
were pooled for a given replicate. 

Altogether, Malpighian tubules were isolated from 1) mosqui- 
toes fed a blood meal at 3 different time points (3 h, 12 h, 24 h; 40 
mosquitoes at each time point) and 2) mosquitoes not fed a blood 
meal at similar time points (40 mosquitoes each). A total of 3 
biological replicates was obtained for each time point in both the 
BF and NBF groups, resulting in 18 sets of tubules for RNA 
isolation and cDNA libraiy preparation (triplicates each of 3 h BF, 
12 h BF, 24 h BF, 3 h NBF, 12 h NBF, and 24 h NBF). 

Isolation and purification of total RNA and assessment of 
its quantity and quality 

Total RNA was extracted from each set of Malpighian tubules 
immediately after they were isolated from the mosquitoes using the 
method of Chomczynski and Sacchi [23]. The resulting RNA was 
treated with TURBO DNA^f« (Life Technologies) to remove 
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genomic DNA and then purified with a RNA Clean & 
Concentrator-5 kit (Zymo Research, Irvine, CA), according to 
the manufacturer's protocol. 

The purified RNA was initially measured for quantity and 
quality using a NanoDrop 2000c Spectrophotometer (Thermo 
Fisher Scientific, Waltham, MA). Samples with a concentration < 
20 ng/|il or poor absorbance ratios (i.e., 260/280 value <1.6; 
260/230 value <1.6 or >3.0) were discarded. RNA quality was 
further assessed using the Experion Automated Electrophoresis 
System (Bio-Rad, Hercules, CA). Only samples with RNA Quality 
indicator values of 7.5 or higher were used. The concentration of 
RNA was determined with a Qubit 2.0 Fluorometer (Life 
Technologies). 

Synthesis and sequencing of cDNA libraries 

Total RNA (565 ng) was used to synthesize adaptor-indexed 
double-stranded cDNA libraries using the TruSeq DNA Sample 
Prep Kit V2, Set A and B (lUumina, San Diego, CA). The size 
chosen for libraries was ~270 bp. The quality of the synthesized 
libraries was evaluated using the Agilent 2100 Bioanalyzer High 
Sensitive DNA Chip (Agilent Technologies, Santa Clara, CA) and 
the quantity determined using the Qubit 2.0 Fluorometer (Life 
Technologies). 

The 18 resulting cDNA libraries were diluted to 18 nM and 
pooled to generate a multiplexed cDNA libraiy (using 18 unique 
indexed adapters) of 36 fM. The pooled library was sequenced 
using the lUumina HiSeq 2000 platform at the Ohio State 
University Comprehensive Cancer Center (Columbus, OH). 
Demultiplexing was performed with CASAVA 1.8.2. FASTQ 
fUes were generated from the 'basecaU' files. All single-end reads 
were submitted to the NCBI sequence read archive (accession 
number SRP034701). The sequencing of all 18 libraries generated 
over 232 million single-end raw reads (~13 million reads per 
sample) (Table SI). 

Filtering data, alignment, and read counting 

The MCIC-Galaxy pipeline was implemented for preprocess- 
ing, filtering, and data analysis [24]. The raw reads were first 
analyzed with the "FASTQC" tool (http://www.bioinformatics. 
babraham.ac.uk/projects/fastqc/) to assess quality. Adapters were 
removed using CUT AD APT [25] with an error rate of 0. 1 and a 
minimum overlap length of 6. Reads were trimmed for length and 
quality using the "Trim the reads" tool, version 1.2.2, (Phred 
threshold score <20; read length <20 bp), discarding all miscalled 
bases, but not duplicates or polyA tails, because reads were aligned 
on a reference transcriptome. The number of reads retained after 
removing adapters, low quality, and short reads was > 1 94 million, 
which is ~84% of the total number of original 232 million reads. 

After the preprocessing and filtering, reads were aligned to the 
following two reference transcriptomes using "Burrow-Wheeler 
Aligner", version 1.2.3 [26]: 1) A. aegypti (www.vectorbase.org; 
AaegLl.4, vl.OO; 18,769 sequences) [27] and 2) A. albopictus (www. 
ncbi.nlm.nih.gov;transcript shotgun archive accession numbers 
J0845359-J09 13491, 68,413 sequences) [28]. The number of 
reads that aligned to unique and redundant transcripts in the 
reference transcriptome was determined for each sample using 
"Cotmt Features" tool, version 0.91. The dataset was filtered to 
contain only transcripts with a minimum of five mapped reads for 
any two replicates in at least one treatment/time combination. 
Only the subset of transcripts meeting this criterion was used in 
subsequent analyses. 

Considerably more reads mapped onto the A. albopictus 
transcriptome (~9.6 million reads/sample; 75% reads mapped) 
compared to the A. aegypti transcriptome (~2.4 million reads/ 



sample; 18.5% reads mapped). However, the A. albopictus 
transcriptome assembly and annotation is incomplete and contains 
redundancies (i.e., there are several transcript identification 
numbers corresponding to the same transcript), which limited 
our ability to accurately identify and quantify specific individual 
transcripts of interest (e.g., aquaporins, NH3 detoxification 
enzymes). Moreover, mapping onto the A. aegfpti transcriptome 
resulted in the detection of 9,8 1 3 non-redundant transcripts, which 
is within the range of ~7, 000-18,000 transcripts detected in RNA- 
seq studies of A. aegypti isolated tissues or whole animals 
[12,17,29,30]. Thus, we use the A. aegypti transcriptome and its 
transcript nomenclature for our downstream analyses. 

For quality assessment, "Count Features" output was imple- 
mented to examine the dispersion of the biological replicates 
(libraries per treatment). Counting reads were normalized 
calculating the RPKM (reads per kilobase per million reads). 
Pearson's correlation and Principal Component Analysis graphics 
were generated to assess the similarity of the replications for each 
condition. 

Transcriptome analysis 

The read counts generated using "Count Features" output were 
submitted to the "DESeq" pipeline, version 1.0.0 [31], to identify 
transcripts with a significant differential expression between 
treatments/time points. This tool is based on the negative 
binomial distribution, with normalized libraries by size factor 
(developing an estimate effective library size) [31]. The compar- 
isons made were a NBF treatment against BF treatments at three 
time points (3 h, 12 h and 24 h). All of the transcripts with a 
significant expression change were filtered by a FDR-adjusted P 
value threshold of 0.05, and their log2 fold-change value was 
recorded. 

The DAVID v6.7 annotation clustering module [32] was used 
to classify difiFerentiaUy expressed transcripts into functional 
groups. Clustering analysis was carried out for the subset of 
transcripts that had showed sustained up- or down-regulation (i.e., 
at least two consecutive time points of differential expression). The 
DAVID is currently not compatible with A. aegypti transcript IDs. 
Thus, the differentially-expressed transcripts were first converted 
to A. gambiae transcript IDs using tBLASTx (E-value <10 ^''). 
Then, enrichment of GO and other annotation terms in candidate 
sub-lists were explored using the functional annotation clustering 
tool. This clustering method condenses the input transcript list into 
functionally related transcripts (annotation clusters), taking into 
account the similarity of their annotation profiles based on 
multiple annotation sources (e.g. GO terms and Interpro key- 
words). 

The clusters are assigned an enrichment score, which represents 
the minus log-transformed geometric mean of the modified Fisher 
Exact (EASE) Scores within the cluster. Significantly enriched 
annotation clusters were defined as those containing a minimum 
enrichment score of 1.3, because the -log (0.05) = 1.3. Thus, an 
enrichment score of >1.3 corresponds to a P<0.05. The 
enrichment score is based on the following parameters: Similarity 
Term Overlap = 3; Similarity Threshold = 0.7; Initial Group 
Membership = 3; Final Group Membership — 5; Multiple 
Linkage Threshold = 0.3. Once the significantly-enriched 
functional clusters were identified, the A. gambiae transcript IDs 
within them were converted back to their respective A. aegypti 
transcript IDs. 

We also performed a preliminary DAVID clustering analysis 
using reads mapped to the existing A. albopictus reference 
transcriptome, despite its incomplete assembly and annotation. 
Notably, there was a good correlation between the enriched 
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functional clusters identified using this reference and those using 
the A. aegypti reference. For example, the thioredoxin, glutathione 
S-transferase, vitamin binding, cofactor metabolic process, oxida- 
tive phosphorylation/ATP synthesis, protein biosynthesis, glycoyl- 
sis, and ATPase activity clusters were enriched in both DAVID 
analyses. Thus, our decision to use the A. aegypti transcriptome, 
which allows for better identification and quantification of specific 
transcripts (see above), does not compromise our ability to identify 
enriched functional pathways. 

Results and Discussion 

As explained in the Methods, Malpighian tubules were isolated 
from mosquitoes at 3 diflFerent time points after being fed a blood 
meal (3 h, 12 h, 24 h). In parallel, Malpighian tubules were 
isolated from mosquitoes that were not fed a blood meal at similar 
time points. A total of 3 biological replicates were obtained for 
each time point in both the blood fed (BF) and non-blood fed 
(NBF) groups. 

A principal component analysis (PC A) of the Malpighian-tubule 
transcriptomes derived from the 18 sequenced cDNA libraries 
(triplicates each of 3 h BF, 12 h BF, 24 h BF, 3 h NBF, 12 h NBF, 
and 24 h NBF) revealed distinct clustering of samples by treatment 
and time (Figure 1). Notably, all of the NBF samples clustered 
tightly with one another, regardless of time. Consistent with this 
observation, a GLM (general linear model) two-way ANOVA 
comparing the RPKM values among the NBF samples revealed no 
significant differences (F-value = 0.731; Table S2), which 
indicates that time does not affect transcript accumulation in the 
NBF tubules. The PCA also showed that the BF samples clustered 
separately from their respective NBF controls and each other 
(Figure 1). Consistent with this observation, an ANOVA 
comparing the RPKM values among the BF samples revealed 
significant diflFereiices among all of the samples (i^-value = 0.017; 
Table S3), which indicates that time affects transcript accumula- 
tion in the BF tubules. 

Below, we dissect the quantitative and qualitative differences 
among the NBF and BF samples in more detail. For simplicity and 
consistency, we selected the 24 h NBF samples as a universal 
control, because 1) an ANOVA found no significant transcript 
differences across all the time points for the NBF samples (see 
above) and 2) the 24 h NBF samples exhibited the lowest internal 
variation among the NBF samples (Figure 1). As described in the 
Methods, we used an assembled transcriptome for A. aegypti as a 
reference given its superior annotation and nominal redundancy 
compared to the available transcriptome for A. albopictus, which is 
stiU in its early stages of development. 

Differential expression in Malpighian tubules after blood 
feeding 

DESeq was used to search for transcripts differentially expressed 
between the 24 h NBF control and the BF treatments at each time 
point. Using the A. aegfpti transcriptome as a reference, a total of 
1,857 non-redundant transcripts was found to be diflferentially 
expressed over all of the time points (~10% of the A. aegypti 
transcriptome). Table 1 shows that the Malpighian tubules from 
the BF mosquitoes were characterized by progressive increases in 
differential expression throughout the time series. At each time 
point, the differentially expressed transcripts consist of similar 
numbers of up- and down-regulated transcripts. 

We next aimed to identify enriched, functional pathways within 
the differentially expressed transcripts using a DAVID functional 
clustering analysis [32,33]. We focused our analysis on transcripts 
that exhibited 'sustained' changes in differential expression after 



blood feeding, which we defined as those significantly up- or 
down-regulated for at least two consecutive time periods. Based on 
the A. aegypti transcriptome, a total of 669 transcripts met our 
'sustained' criterion, consisting of 340 up-regulated transcripts and 
329 down-regulated transcripts. As shown in Table 2, the DAVID 
analysis revealed a significant enrichment (enrichment score > 
1.3) of 1) nine functional groups among the sustained, up-regulated 
transcripts and 2) six functional groups among the sustained, 
down-regulated transcripts. The identities of the transcripts that 
comprise the functional groups of Table 2 and their respective 
heat maps of differential expression are shown in Figures S1-S15. 

Cursory interpretations of the changes to these broad functional 
groups and the transcripts within them suggests that blood feeding 
promotes the expression of transcripts associated with 1) antiox- 
idants and detoxification, 2) proteolytic activity, 3) amino acid 
metabolism, and 4) cytoskeletal dynamics. On the other hand, 
blood feeding appears to suppress the expression of transcripts 
associated with 1) oxidative metabolism, 2) active transport, and 3) 
mRNA translation. Below, we discuss these interpretations in 
more detail with the caveat that transcript levels may not 
necessarily reflect protein abundance, biochemical activity, or 
physiological function. Thus, we consider our interpretations as 
the buUdiiig of hypotheses that will require testing in future studies 
using functional genetic, biochemical, and physiological tech- 
niques. 

Down-regulation of transcripts associated with oxidative 
metabolism and active transport — Potential implications 
for diuretic fluid secretion 

The mosquito Malpighian tubule epithelium is well-studied, 
because of its remarkable capacity for active transepithelial fluid 
secretion, which mediates the post-prandial diuresis. The vacuolar 
(V-type) H"*"- ATPase is the ultimate energizer of transepithelial 
fluid secretion in the epithelium [34] . This proton pump resides in 
the luminal brush border of principal cells where it is situated in 
close proximity to mitochondria that fuel the pump with ATP 
[35,36]; the pump is a multisubunit protein consisting of two 
sectors: 1) a catalytic, cytosolic Vi sector and 2) a H'^-translocating, 
membrane-bound Vq sector [37]. Inhibiting the pump or the 
production of ATP in the epithelium effectively inhibits fluid 
secretion [34,38]. Thus, the enrichment of the 'oxidative 
phosphorylation/ATP synthesis', 'ATPase activity', 'glycolysis', 
and 'sugar/inositol transporter' functional clusters among the 
transcripts that exhibited a sustained down-regulation (Table 2) 
caught our attention. 

Listed prominently among the down-regulated transcripts in the 
'oxidative phosphorylation/ ATP synthesis' (Figure SI) and 
'ATPase activity' (Figure S2) functional clusters are those encoding 
subunits of the V-type H"'"-ATPase. Remarkably, as shown in 
Figure 2, fourteen transcripts encoding subunits of the V-type H"*"- 
ATPase exhibit a sustained down-regulation after blood feeding, 
while only one shows a sustained up-regulation (AAEL003743- 
RA). Furthermore, a manual search of all the differentially- 
expressed transcripts (including those not considered 'sustained') 
revealed three other transcripts associated with the V-type FF^- 
ATPase that are down-regulated transiently at one or two non- 
consecutive time points (Le., AAEL002464-RA, AAEL010819- 
RA, AAEL010819-RB in Figure 2). Among all of these down- 
regulated transcripts, nine encode subunits of the Vi sector of the 
V-type H^-ATPase, seven encode subunits of the Vo sector, and 
one encodes an accessory protein (Figure 2). The only up- 
regulated transcript encodes subunit 'a' of the Vq sector (a.k.a. 
vhalOO-1). 
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Figure 1. Principal component analysis of transcript expression in Malpigliian tubules isolated from non-blood fed mosquitoes 
(NBF) at 3 h (gray circles), 12 h (white circles), and 24 h (black circles) and blood-fed (BP) mosquitoes at 3 h (blue circles), 12 h (red 
circles), and 24 h (green circles). Principal component 1 is time and principal component 2 is treatment (NBF, BF). 
doi:1 0.1 371 /journal.pntd.0002929.g001 



The above changes to transcripts of the V-type H'^'-ATPase in 
the Malpighian tubules of A. albopictus contrast with those 
previously reported in the midgut of ^. aegypti after blood feeding 
[19]. In our study, nearly aU subunits of the V-type H^-ATPase 
exhibited a down-regulation in Malpighian tubules after blood 
feeding (Figure 2), whereas Sanders et al. found that transcripts 
encoding V-type H"'"-ATPase subunits exhibited an up-regulation 
in the midgut of ^. aeg)!pti at 12 h and 24 h after blood feeding 
[19]. These contrasting results suggest that the regulation of V- 
type H"''-ATPase expression in response to blood feeding is tissue 
dependent in mosquitoes. 

Notable among the transcripts that exhibited a sustained down- 
regulation after blood feeding in the 'oxidative phosphorylation/ 
ATP synthesis' (Figure SI) and 'glycolysis' (Figure S3) functional 
clusters are those encoding enzymes associated with glycolysis, the 
citric acid cycle, and ATP synthesis, such as phosphofructokinase, 
pyruvate kinase, enolase, 2-oxoglutarate dehydrogenase, glycerol- 
3-phosphate dehydrogenase, and acetyl-CoA synthetase. In 
addition, within the 'oxidative phosphorylation/ ATP synthesis' 
(Figure SI) and 'sugar/inositol transporter' (Figure S4) functional 
clusters are transcripts encoding putative SLC2-like sugar trans- 
porters, which import glucose into cells for use as a fuel to generate 
ATP. Thus, in Malpighian tubules, blood-feeding leads to a 
decrease in the abundance of transcripts encoding enzymes 



associated with the catabolism of glucose and synthesis of ATP, 
which is consistent with the aforementioned decrease in abim- 
dance of transcripts encoding subunits of the V-type H'*"-ATPase 
subunits. 

Also notable among the transcripts listed in the 'oxidative 
phosphorylation/ ATP synthesis' (Figure SI) functional cluster are 
those encoding ion transport mechanisms that are known or 
hypothesized to play a role in the transepithelial secretion of ions by 
mosquito Malpighian tubules. We discuss these mechanisms below. 

Mechanisms of uptake. Both barium-sensitive inward- 
rectifying (Kir) channels and bumetanide-sensitive Na,K,2Cl- 
cotransporters (NKCCs) are considered the major mechanisms for 
the uptake of K'*' from the hemolymph across the basolateral 
membranes of mosquito Malpighian tubules [39,40,41], where 
they presumably localize to principal cells. As shown in Figure 3, 
blood-feeding ehcits a sustained down-regulation of two transcripts 
encoding Kir channel subunits (AAEL008931-RA, AAELO 13373- 
RA) and one transcript encoding a putative NKCC 
(AAEL006180-RA). We recently cloned and characterized 
AAEL008931-RA (Kir2B) from the Malpighian tubules of A. 
aegypti and have confirmed that it encodes a barium-sensitive Kir 
channel [41]. The NKCC is an ortholog oi Drosophila melanogaster 
Ncc69, which has been shown to a play a key role in fluid and K""" 
secretion in fruit fly Malpighian tubules [42] . These changes in Kir 
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Table 1. Total number of transcripts whose abundance Is 
significantly affected by blood feeding at each time point 
compared to non-blood fed controls. 



Treatment 


3 h BF 


12 h BF 


24 h BF 


Dp-regulation 


207 


553 


564 


Down-regulation 


192 


482 


669 


Total 


399 


1035 


1233 



doi:10.1371/journal.pntd.0002929.t001 



and NKCC transcript abundance suggest that blood-feeding leads 
to a decrease in the capacity of the tubules to uptake K*^ across the 
basolateral membrane, which is expected to reduce the diuretic 
capacity of the tubules. 

Mechanisms of Na^ uptake. A cAMP-induced, amUoride- 
sensitive Na"''-conductance [40,43,44] and an amUoride-insensitive 
Na/H exchanger (NHE3) [45] are considered important mecha- 
nisms for the uptake of Na""" across the basolateral membranes of 
mosquito Malpighian tubules, where they are hypothesized and 
known, respectively, to localize to principal cells. Although the 
molecular identity of the Na''"-conductance is unknown, one 
transcript encoding a putative 'pickpocket' Na""" channel exhibits a 
sustained down-regulation after blood feeding (Figure 3). Pickpocket 
genes contribute to 1) transepithelial movements of fluid in the 
tracheal tubes [46] and 2) osmosensation in the nervous system of 
Drosophila [47]. Thus, it is reasonable to hypothesize a potential role 
of pickpocket in fluid secretion by mosquito Malpighian tubules. 
Moreover, a transcript encoding NHE3 shows a similar pattern of 
down-regulation as the pickpocket transcript after blood feeding 
(Figure 3). These changes in pickpocket and NHE3 abundance 
suggest that blood-feeding leads to a decrease in the capacity of the 
tubules to uptake Na across the basolateral membrane, which is 
expected to reduce the diuretic capacity of the tubules. 

Mechanisms of HCO3" transport. The basolateral mem- 
branes of mosquito Malpighian tubules possess a DIDS-sensitive 
Cl/HCO-j anion exchanger (AE), which localizes to stellate cells 



[48,49]. We have shown that this AE plays a key role in diuretic 
fluid secretion in isolated Malpighian tubules of A. aegppti by 
presumably regulating the oxidative metabolism of principal cells 
[7,48]. The AE may also contribute to the transepithelial secretion 
of Cr by stellate cells [7] . Notably, after blood feeding, a transcript 
encoding the AE (AAEL001269-RA) is down-regulated (Figure 3). 
This change suggests that blood feeding leads to a decrease in the 
capacity of stellate cells to regulate the oxidative metabolism of 
principal cells and to transport CI , which is expected to decrease 
the diuretic capacity of the tubules. 

Mechanisms of water transport. Given the changes to the 
above ion transport mechanisms, we also manually searched for 
diflerentially-expressed transcripts encoding aquaporin water 
channels (AQPs). Although AQP transcripts were not part of the 
enriched functional clusters identified by the DAVID analysis, the 
Malpighian tubules of A. aegfpti express at least three AQPs 
(AQPl, AQP4, AQP5) that play an important role in the post- 
prandial diuresis and exhibit differential expression after blood 
feeding [22]. Furthermore, AQPl immunoreactivity is expressed 
in stellate cells and some principal cells in the Malpighian tubules 
of A. gambiae [50,51]. 

As shown in Figure 3, transcripts encoding AQPl and AQP4 
exhibit a sustained down-regulation after blood feeding and a 
transcript encoding AQP2 is down-regulated by 24 h. In contrast, 
two transcripts encoding AQP5 isoforms are up-regulated, but only 
transientiy at 12 h (Figure 3). Surprisingly, the changes in AQP 
expression we observed in A. albopictus Malpighian tubules contrast 
with those found in the Malpighian tubules of A. aegppti by Drake 
and colleagues [22]. In particular, at 3 h after a blood meal, AQPl, 
AQP2, and AQP4 each exhibited transient increases of expression 
in the Malpighian tubules of A. aegfpti [22], whereas we observed 
decreases of all three within 24 h in ^. albopictus. Moreover, in A. 
albopictus, we observed a transient increase of AQP5 at 12 h 
(Figure 3), whereas in A. aegypti, AQP5 expression was increased in 
the tubules at 3 h, 12 h, and 24 h after blood feeding [22]. These 
contrasts may point to species-specitic differences in how the AQPs 
of mosquito Malpighian tubules respond to the physiological stresses 
encountered during blood feeding. 

The consensus of the above findings in Figures S1-S4 and 
Figures 2-3 suggests that as early as 3 h after a blood meal, the 



Table 2. DAVID functional clusters that are significantly enriched among the transcripts in Malpighian tubules exhibiting a 
sustained, up-regulation or down-regulation after blood feeding. 





Functional cluster 


Number of 
Transcripts 


Enrichment 
Score 


Functional cluster 


Number of 
Transcripts 


Enrichment 
Score 


Up-regulated 






Down-regulated 






Thioredoxin 


17 


4.05 


Oxidative phosphorylation/ATP 
synthesis 


40 


5.28 


Cofactor metabolic process 


16 


1.76 


Protein biosynthesis 


23 


6.19 


Tubulin, GTPase 
domain 


15 


1.61 


Glycolysis 


22 


5.52 


Amine biosynthetic process 


11 


1.86 


ATPase activity 


13 


1.83 


ATPase, AAA+ type 


9 


1.42 


Translatlonal elongation 


6 


3.07 


Glutathione 
S-transferase 


8 


1.55 


Sugar/lnositol transporter 


5 


2.14 


Vitamin binding 


8 


1.7 








Vitamin biosynthetic process 


8 


13.1 








Proteasome complex 


7 


3.05 









The transcripts that comprise each functional group are listed In Figures S1-S15. 
doi:1 0.1 371 /journal.pntd.0002929.t002 
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Transcript ID 


Sector, subunit 


AAEL008787-RA 


Vi, A 

1 ' 


AAEL005798-RA 


Vi, B 


AAEL005173-RA 


Vi, c 


AAEL015594-RA 


Vi, c 


AAEL009808-RA 


Vi, D 


AAEL012035-RA 


Vi, E 


AAEL002464-RA 


Vi, F 

1 ' 


AAEL012819-RA 


Vi, G 

1 ) 


AAEL006516-RA 


Vi, H 

1 ' 


AAEL006390-RA 


Vo, a 


AAEL014053-RA 


Vo, a 


AAEL003743-RA 


Vo, a 


AAEL000291-RA 


Vo, c 


AAEL012113-RA 


Vo, c" 


AAEL011025-RA 


Vo, d 


AAEL010819-RA 


Vo,e 


AAEL010819-RB 


Vo,e 


AAEL007777-RA 


Accessory, S1 



3h 


12 h 


24 h 




































1 












-5 


0 5 



Figure 2. Differential expression after blood feeding of transcripts encoding subunits of the V-type H^-ATPase. Blue shading indicates 
significant down-regulation compared to NBF controls, whereas yellow shading indicates significant up-regulation compared to NBF controls. Lack of 
shading indicates no significant difference relative to NBF controls. Degree of shading is based on value of Log2 fold change as indicated by the scale 
below on the right. 
doi:1 0.1 371 /journal.pntd.0002929.g002 



Malpighian tubule epithelium initiates a host of molecular changes 
that are expected to reduce its capacity for 1) generating 
transepithelial electrochemical gradients by the V-type H^-ATPase, 
and 2) mediating the transepithelial transport of K"*", Na"*", CI and 
water. Thus, during the chronic processing of blood meals (24—48 h 
post blood meal) the capacity of the Malpighian tubules to mediate 
active, transepithelial fluid secretion (i.e., diuresis) may decrease. 
Such a reduction in diuretic capacity would make physiological 
sense given that the major stresses to hemolymph salt and water 
homeostasis are alleviated within 2 h after blood feeding [8] . 

Down-regulation of transcripts associated with the 
translation of mRNA 

The other functional clusters enriched among the transcripts 
that exhibited a sustained down-regulation are related to the 
translation of mRNA (i.e., protein biosynthesis and translational 
elongation; Table 2). These transcripts consist primarily of 
ribosomal protein subunits and translation/ elongation factors that 



exhibit a down-regulation at 12 h and 24 h after a blood meal 
(Figure S5 and Figure S6), which suggests that most newly- 
translated proteins in the tubules in response to blood feeding are 
synthesized within 12 h. Furthermore, the data suggest that during 
the chronic processing of blood meals (24—48 h after blood feeding) 
the capacity of the tubules for de novo protein synthesis may 
decrease. These observations in the Malpighian tubules of A. 
alhopictus are similar to results of pre\'ious studies in the midgut and 
fat body of ^. aegypti, where blood feeding led to a down-regulation 
of transcripts encoding ribosomal protein subunits and/or transla- 
tion factors within 12 h to 24 h after a blood meal [17,19,52]. 

Up-regulation of transcripts associated with antioxidants 
and detoxification — Implications for heme metabolism/ 
detoxification and metabolite excretion 

Hemoglobin is the most abundant protein in mammalian blood, 
and its digestion leads to the production of heme, a highly toxic 
metabolite that causes cell and tissue damage via oxidative stress 
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Transcript ID 


Predicted encoded protein (abbreviation) 


AAEL008931-RA 


IMWdlU-ltJuUIVIilLj r\ -L/lldilllcl IrxllZD) 


AAEL013373-RA 


Kir2B' 


AAEL006180-RA 


Na.K.CI cotransporters (NKCC1) 


AAEL006613-RA 


pickpocket sodium channel (ppk) 


AAEL001503-RA 


Na,H exchanger (NHE3) 


AAEL001269-RA 


Anion exchanger (AE) 


AAEL003512-RA 


Aquaporin (AQP1) 


AAEL003550-RA 


AQP2 


AAEL005001-RA 


AQP4 


AAEL005008-RA 


AQP5 


AAEL005008-RC 


AQP5 




-5 



0 



Figure 3. Differential expression after blood feeding of transcripts encoding known and putative ion and water transport 
mechanisms. Blue shading indicates significant down-regulation compared to NBF controls, whereas yellow shading indicates significant up- 
regulation compared to NBF controls. Lack of shading indicates no significant difference relative to NBF controls. Degree of shading is based on value 
of Log2 fold change as indicated on the right. 
doi:1 0.1 371/journal.pntd.0002929.g003 



and/or the disruption of plasma membranes [53]. Mosquitoes 
utilize a variety of mechanisms to detoxify heme and limit its 
absorption into the hemolymph. The first line of defense is found 
in the midgut, which secretes a peritrophic matrix that encapsu- 
lates the ingested blood cells and may sequester more than half of 
the amount of heme in a typical blood meal [54,55] . Additional 
mechanisms for heme detoxification include the I) enzymatic 
degradation of heme by heme oxygenase (HO) to produce 
bHiverdin, Fe^^, and carbon monoxide [56], 2) chelation of heme 
by xanthurenic acid PCA), which is a product of the kynurenine 
pathway of tryptophan catabolism [57,58], and 3) binding and/or 
catabolism of heme by glutathione S-transferases (GSTs) [59,60] . 
Furthermore, protection against heme-induced, free-radical dam- 
age in mosquitoes can be mediated by: 1) antioxidant enzymes, 
such as glutathione peroxidase (GP), thioredoxin peroxidase 
(THP), thioredoxin reductase (THR), superoxide dismutase 
(SOD), and catalase (CAT); 2) antioxidant proteins, such as 
thioredoxin (TH); and 3) small molecule antioxidants, such as 
glutathione and uric acid [53]. 

Thus, we were intrigued by the enrichment of the 'thioredoxin' 
and 'glutathione S-transferase' functional clusters among the 
transcripts that exhibited a sustained up-regulation after blood 
feeding (Table 2; Figures S7-S8). Furthermore, we noticed that 
the 'cofactor metabolic process' (Figure S9) and 'vitamin 
biosynthetic process' (Figure SIO) functional clusters contained 
several transcripts associated with putative antioxidant and 
detoxification mechanisms, and that the 'ATPase/AAA-F type' 
(Figure S 1 1 ) functional cluster contained several transcripts 
encoding putative ATP-binding cassette (ABC) transporters, which 
play key roles in insect metabolite /xenobiotic excretion [61]. 
Below, we discuss the up-regulation of transcripts after a blood 
meal within the aforementioned functional clusters in the context 
of 1) the prevention of heme-induced oxidative cell and tissue 
damage and 2) the detoxification/ excretion of heme and heme- 
related metabolites. 

Heme detoxification and antioxidant production 
mechanisms. As shown in Figure 4, we fmd a remarkable. 



sustained up-regulation after blood feeding of several transcripts 
encoding putative heme-detoxilication mechanisms, including 4 
transcripts that encode enzymes associated with XA production 
(kynurenine 3-monoxygenases and a tryptophan dioxygenase) and 
8 transcripts encoding GSTs (Delta, Sigma, EpsHon, and Rho 
classes). A manual search for other GSTs among the transcripts 
that were diflFerentiaUy expressed (non-sustained) revealed 3 others: 
1) a putative microsomal GST (AAEL006829-RA), which was 
up-regulated, albeit transiently, at 3 h and 24 h after blood 
feeding, 2) the heme-binding GSTX2 (AAEL010500-RA), which 
was up-regulated transiently at 3 h, and 3) a Zeta class GST 
(AAEL011934-RA), which was down-regulated by 24 h after 
blood feeding (Figure 4). Interestingly, two of the up-regulated 
GST transcripts (AAEL007962-RA, GSTe4; AAEL010500-RA, 
GSTX2) were also found to be up-regulated in pyrethroid- 
resistant hues of ^. aegypti [62,63], which suggests putative roles of 
these GSTs in both heme and insecticide detoxification in 
mosquitoes. 

Transcripts encoding putative antioxidant mechanisms for 
protection against heme-induced oxidative stress are also in general 
up-regulated in Malpighian tubules after blood feeding. Namely, 
transcripts encoding CAT, GP, TH, THP, and THR exhibited a 
sustained up-regulation after a blood meal (Figure 4). A manual 
search for additional differentially-expressed transcripts (non- 
sustained) encoding THR revealed two others (AAEL002886- 
RA,-RB), which were transiendy up-regulated at 3 h (Figure 4); a 
similar search for transcripts encoding THPs found one 
(AAEL004112-RA) that is up-regulated by 24 h and another 
(AAEL013528-RA) that is down-regulated by 24 h (Figure 4). 
Similar patterns in the expression of CAT and THR have been 
observed in the midgut of ^. aeg)ipli after a blood meal [19]. 

The consensus of the above results suggests that the molecular 
capacity of the Malpighian tubules to detoxify heme and protect 
against heme-induced free radical damage increases after mosqui- 
toes ingest a blood meal. Thus, the data may point towards a 
putative new role of the tubule epithelium after a blood meal in the 
detoxification of heme that enters the hemolymph. 



PLOS Neglected Tropical Diseases | www.plosntds.org 



8 



June 2014 | Volume 8 | Issue 6 | e2929 



Changes to the Renal Transcriptome of Mosquitoes after a Blood Meal 



Transcript ID 


Predicted encoded protein (abbreviation) 


AAEL013407-RA 


Catalase 


AAEL012069-RB 


Glutathione peroxidase 


AAEL001061-RA 


Glutatliione transferase (GSTD1 ) 


AAEL001061-RB 


Glutathione transferase (GSTD1 ) 


AAEL001061-RC 


Glutathione transferase (GSTD1 ) 


AAEL007962-RA 


Glutathione transferase (GSTe4) 


AAEL011741-RA 


Glutathione transferase (GSTS1) 


AAEL011741-RB 


Glutathione transferase (GSTS1) 


AAEL009017-RA 


Glutathione transferase (GSTt1 ) 


AAEL010500-RA 


Glutathione transferase (GSTX2) 


AAEL011934-RA 


Glutathione transferase (GSTzl) 


AAEL006764-RA 


Glutathione transferase D10, putative (GSTD10) 


AAEL008879-RA 


Kynurenlne 3-monooxygenase 


AAEL008879-RB 


Kynurenlne 3-monooxygenase 


AAEL012757-RA 


Kynurenlne 3-monooxygenase 


AAEL000428-RA 


Tryptophan 2,3-dioxygenase 


AAEL006829-RA 


Microsomal glutathione s-transferase 


AAEL010777-RA 


Thioredoxin (TRX), putative 


AAEL002309-RA 


Thioredoxin peroxidase 


AAEL004112-RA 


Thioredoxin peroxidase 


AAEL013528-RA 


Thioredoxin peroxidase 


AAEL014548-RA 


Thioredoxin peroxidase 


AAEL002886-RA 


Thioredoxin reductase 


AAEL002886-RB 


Thioredoxin reductase 



3 h 



12 h 



24 h 



-5 0 5 

Figure 4. Differential expression after blood feeding of transcripts associated with heme detoxification and antioxidant 
mechanisms. Blue shading indicates significant down-regulation compared to NBF controls, whereas yellow shading indicates significant up- 
regulation compared to NBF controls. Lack of shading indicates no significant difference relative to NBF controls. Degree of shading is based on value 
of Log2 fold change as indicated by the scale below on the right. 
doi:1 0.1 371 /journal.pntd.0002929.g004 



Heme and heme-metabolite excretion. In addition to 
direct heme detoxification, tlie up-regulation of GSTs can be 
interpreted as serving a potential role in improving the clearance 
of heme metabolites, such as biliverdin (derived from HO activity), 
from the hemolymph. Mammalian GSTs possess binding-sites for 
various organic molecules, including biliverdin, bilirubin, and even 
heme [64,65,66]. In the mammalian Uver, the bincUng of bilirubin 
by intracellular GSTs provides a favorable diffusional gradient for 
bilirubin uptake by liver cells [65] . Moreover, GSTs are known to 
conjugate reduced glutathione onto large organic molecules 
thereby increasing their solubility, decreasing their toxicity, and 
enhancing their transport by ABC transporters [67]. Thus, it is 
reasonable to hypothesize that one or more of the up-regulated 
GSTs in Malpighian tubules may contribute towards 1) maintain- 
ing a favorable chemical gradient for the uptake of heme and/ or 
biliverdin from the hemolymph, and/or 2) chemically modifying 
heme and/ or biliverdin to improve their excretion. 

In regards to the potential excretion of heme and heme-derived 
metabolites by Malpighian tubules, putative ABC transporters are 
among those in the 'ATPase/AAA+ type' functional cluster 
(Figure Sll). In general, ABC transporters are promiscuous 
in their transport substrates, which include heme, bilirubin, 



glutathione-conjugated organic compounds, uric acid, xenobiotics 
(including insecticides), ions (e.g., cytstic-fibrosis transmembrane 
regulator), and lipids [68,69,70,71,72]. As shown in Figure 5, five 
transcripts encoding ABC transporters (of the 'A', 'B', and 'C 
types) exhibit a sustained up-regulation after blood feecUng. 
Furthermore, a manual search of all the differentially-expressed 
transcripts (non-sustained) revealed eight others (of the 'A', 'B', 'C, 
and 'G' types) that exhibited a transient up-regulation at one or 
two discontinuous time points (Figure 5). Only one transcript 
(AAEL005929-RA) exhibited a down-regnlation and another 
(AAEL005937-RA) exhibited an up-regulation at 3 h followed 
by a down-regulation at 12 h. Thus, there is a general up- 
regulation of transcripts encoding ABC transporters, which 
suggests that blood feeding leads to an increased capacity of the 
tubule epithelium for the excretion of metabolites and xenobiotics 
derived from the blood meal. 

Up-regulation of transcripts associated with 
proteasomes — Implications for proteolytic activity 

Another functional cluster enriched among the transcripts that 
exhibited a sustained up-regulation is the 'proteasome complex' 
(Table 2), which is related to the degradation of protein. These 
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Transcript ID 


Predicted encoded protein (abbreviation) 


AAEL014d99-RA 


ABC Transporter family 'A' 


A Am OHOcriO da 


ABC Transporter family 'A' 


AA^i r\r\Q A A r^A 


ABC Transporter family 'B' 


AA^i r\ A r\'~i~7 r\ nA 

AAEL01 0379-KA 


ABC Transporter family B 


A Am r\r\ A~7 A n da 

AALL(J(J474o-KA 


ABC 1 ransporter family C 


AA^I r\r\ A~7 A'^ l~) [-) 

AAEL00474O-KB 


ABC Transporter family 'C 


AAcri ^\^o^o^ da 
AAbLUIzl oy-KA 


ABC 1 ransporter family C 


AA^i r\ncoon da 


ABC Transporter family 'C 


AAELUUbyo / -KA 


ABC Transporter family 'C 


AA^I AHOOOO OA 

AAELUIzooD-KA 


ABC 1 ransporter family C 




ABC Transporter family 'C 


AAEL013854-RA 


ABC Transporter family 'C 


AAEL008629-RA 


ABC Transporter family 'G' 


AAEL008631-RA 


ABC Transporter family 'G' 


AAEL008635-RA 


ABC Transporter family 'G' 



3h 



12 h 



24 h 



Figure 5. Differential expression after blood feeding of transcripts encoding putative ABC transporters. Blue shading indicates 
significant down-regulation compared to NBF controls, whereas yellow shading indicates significant up-regulation compared to NBF controls. Lack of 
shading indicates no significant difference relative to NBF controls. Degree of shading is based on value of Log2 fold change as indicated by the scale 
below on the right. 
doi:1 0.1 371 /journal.pntd.0002929.g005 



transcripts consist entirely of those encoding proteasome or 
protease regulatory subunits and most exhibit an up-regulation 
at 3 h and 12 h after a blood meal (Figure SI 2), which suggests 
that the tubules enhance their molecular capacity for proteolytic 
activity early after a blood meal. Enhanced proteolytic activity 
within the tubule epithelium may lead to the degradation of 1) 
proteins encoded by the transcripts that are down-regulated after 
blood feeding, and/or 2) proteins that may experience oxidative 
damage from heme. An increase of proteasome activity would also 
be expected to result in an increase of free amino-acids available 
for the synthesis of new proteins [73] — ^perhaps those associated 
with the 'thioredoxin' and 'glutathione S-transferase' functional 
clusters mentioned above. 

Up-regulation of transcripts associated with amino acid 
metabolism — Potential insights into ammonia 
detoxification 

Consistent with the enrichment of the 'proteasome complex' 
functional cluster among up-regulated transcripts, which is 
expected to increase the abundance of free amino acids (see 
above), there is a corresponding enrichment in the 'amine 
biosynthetic process' cluster (Table 2). The transcripts within this 
cluster consist primarily of those encoding enzymes associated with 
amino-acid catabolism and/ or biosynthesis, such as cysteine 
dioxygenase, 2-amino-3-ketobutyrate coenzyme A ligase, gluta- 
mine synthetase, phosphoserine phosphatase, ornithine decarbox- 
ylase, and phosphoserine aminotransferase (Figure SI 3). Similar 
(as well as redundant) transcripts are also found in the 'vitamin 
binding' functional cluster, such as alanine-glyoxylate aminotrans- 
ferase and alanine aminotransferase (Figure SI 4). Notably, the up- 
regulation of these transcripts occurs at 12 li and 24 h, appearing 
to foUow the up-regulation of transcripts associated with the 



proteasome (Figure SI 2). Thus, the potential increase in the 
availability of amino acids derived from proteasome activity in the 
tubules may be followed by an increased molecular capacity to 
breakdown and/or convert amino acids into other products. 

The aforementioned up-regulation of glutamine synthetase and 
alanine aminotransferase drew our attention to a potential role of 
the Malpighian tubules in the handling of ammonia. As 
mosquitoes metabolize a protein-rich blood meal, they face a 
potentially toxic accumulation of ammonia in their hemolymph 
and tissues from the catabolism of proteins and amino acids in the 
blood meal. Glutamine synthetase and alanine aminotransferase 
play prominent roles in detoxifying the ammonia (see below and 
Figure 6). 

In brief, to prevent the build-up of toxic ammonia, mosquitoes 
have a remarkable capacity to fix and assimilate it into free amino 
acids (e.g., alanine, proline, glutamine) via a series of biochemical 
reactions catalyzed by enzymes, such as glutamine synthetase (GS), 
glutamate dehydrogenase (GDH), glutamate synthase (GltS), 
alanine aminotransferase (ALAT), and pyrrolidine-5-carboxylate 
synthase (P5CS) and reductase (P5CR) (Figure 6) [9,74,75,76]. 
Moreover, glutamine can serve as a substrate for the production of 
uric acid through a pathway that includes xanthine dehydrogenase 
(XDH) among other enzymes [9,76]. Once uric acid is produced, 
it can be excreted directly or further converted into allantoin, 
allantoic acid, and urea through a series of biochemical reactions 
catalyzed by urate oxidase (UO), allantoinase (ALLN), and 
aUantoicase (AALC), respectively (Figure 6) [9,76]. 

Figure 7 shows the transcripts associated with this pathway in 
Malpighian tubules that are differentially expressed after blood 
feeding. Namely, two transcripts encoding GS, one transcript 
encoding ALAT, and one transcript encoding XDH exhibit a 
sustained up-regulation after blood feeding. Three transcripts 
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Figure 6. Biochemical pathways for ammonia detoxification in mosquitoes. Redrawn from Isoe and Scaraffia [9], Yellow and blue text 
indicates an enzyme encoded by a transcript that is up-regulated and down-regulated, respectively, after blood feeding. 
doi:1 0.1 371 /journal.pntd.0002929.g006 



encoding a GDH and two transcripts encoding ALAT are each 
transiently up-regulated at 1 2 h, whereas one transcript encoding 
GltS is transiently down-regulated at 12 h (Figure 7). The up- 
regulation of these transcripts occurs at 1 2 h and/ or 24 h after 
blood feeding, which coincides with the putative availability of 
amino acids from increased proteasome activity in the tubule (see 
above), as well as a period of intense protein digestion of the blood 
meal in the midgut [10, 1 1]. 

Interpreting these molecular findings in the context of the 
ammonia detoxification pathway suggests that the Malpigliian 
tubules detoxify ammonia by converting it to glutamate via GDH 
and/or glutamine via GS (Figure 6). Given the up-regulation of 
ALAT, it is reasonable to propose that any newly-formed 
glutamate is converted into alanine (Figure 6). This putative 
handling of ammonia is similar to that reported for the midgut of 
A. aegypti, which fixes and assimilates ammonia into glutamine 
and alanine, as opposed to the fat body which fixes and 
assimilates ammonia into glutamine and proline [7 7] . The fate of 
the alanine in Malpigliian tubules is unknown, but it is possible 
that it is shuttled to the fat body for conversion into proline, 
which can then be shuttled to the flight muscle for use as an 
energy source [78]. 

The fate of glutamine in Malpighian tubules is also 
unknown, but it is possible that it is converted into uric acid 
for excretion, as indicated by 1) the sustained up-regulation of 
XDH, and 2) the transient down-regulation of GltS, after blood 
feeding (Figure 7). Since transcripts encoding UO, ALLN, and 
ALLC were not differentially expressed after blood feeding 



(data not shown), the uric acid may be directly secreted by the 
tubules for excretion (perhaps by an ABC transporter in 
Figure 5), or it may be retained by the epithelium for use as an 
antioxidant to combat potential oxidative damage due to heme. 
Consistent with the former notion, uric acid is excreted by 
mosquitoes after a blood meal [79]. Likewise, uric acid is 
excreted by tsetse flies and reduviid bugs following a blood 
meal [80,81,82]. 

Up-regulation of transcripts associated with cytoskeletal 
dynamics 

The remaining functional cluster enriched among the tran- 
scripts that exhibited a sustained up-regulation is related to 
cytoskeletal dynamics (i.e., 'tubulin, GTPase domain') (Table 2). 
The transcripts in this functional cluster consist primarily of those 
encocUng components of microtubules, such as tubulin chains 
(alpha and beta), dynein light chains, and microtubule-associated 
proteins (Figure SI 5). Another related transcript populating this 
cluster is one encoding a putative Rab GTPase (AAEL006091- 
RA); Rab GTPases play important roles in regulating microtu- 
bule-mediated trafficking of vesicular cargo [83]. These changes 
suggest that blood feeding leads to a more dynamic niicrotubule- 
based cytoskeleton, which may be associated with the intracellular 
trafficking of vesicles and/ or organelles. At least one study has 
found that the actin cytoskeleton of principal cells plays a key role 
in modulating diuretic fluid secretion in A. aegppti Malpighian 
tubules [84]. 
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Transcript ID 


Predicted encoded protein (abbreviation) 


AAEL001887-RA 


Glutamine synthetase (GS) 


AAEL001887-RB 


GS 


AAEL014768-RA 


Glutamate synthase (GItS) 


AAEL010464-RA 


Glutamate dehydrogenase (GDH) 


AAEL010464-RB 


GDH 


AAEL010464-RC 


GDH 


AAEL009872-RA 


Alanine aminotransferase (ALAT) 


AAEL009875-RA 


ALAT 


AAEL009875-RB 


ALAT 


AAEL010630-RA 


Xanthine dehydrogenase (XDH) 



3h 



12 h 



24 h 



-5 0 5 

Figure 7. Differential expression after blood feeding of transcripts associated with ammonia detoxification. Blue shading indicates 
significant down-regulation compared to NBF controls, whereas yellow shading indicates significant up-regulation compared to NBF controls. Lack of 
shading indicates no significant difference relative to NBF controls. Degree of shading is based on value of Log2 fold change as indicated by the scale 
below on the right. 
doi:1 0.1 371/journal.pntd.0002929.g007 
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Figure 8. Summary of the molecular changes identified in the Malpighian tubule epithelium after blood feeding in the present 
study, with an emphasis on the mechanisms of transepithelial fluid secretion. Blue-shaded items represent putative mechanisms or 
pathways that are down-regulated after blood feeding, whereas yellow-shaded text boxes represent putative pathways or functions that are up- 
regulated after blood feeding. See text for details. Gray-shaded items are ion transport mechanisms not examined in the present study. AE, anion 
exchanger [48]; AQP, aquaporin [22]; ATP, adenosine triphosphate [38]; ADP adenosine diphosphate; CIC, chloride channel [86]; Kir, inward-rectifying 
potassium channel [41]; KCC, potassium-chloride cotransporter [87]; PPK, pickpocket sodium channel (hypothesized, present study); NHA, sodium- 
hydrogen antiporter [88]; NHE, sodium-hydrogen exchanger [45]; NKCC, sodium-potassium-chloride cotransporter [89]; VHA, V-type H*-ATPase [35]. 
doi:1 0.1 371 /journal.pntd.0002929.g008 
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Putative changes in tlie microtubule-based cytoskeleton would 
also be consistent with a potential functional transition of the 
epithelium after blood feeding. For example, if the capacity of 
the tubule for diuresis indeed decreases, as expected by the 
down-regulation of transcripts associated with the V-type H"**- 
ATPase and other ion/water transport mechanisms (Figures 2-3), 
then enhanced vesicular trafficking may play a key role in the 
endocytosis and degradation of these membrane-bound proteins. 
Likewise, if the capacity of the tubules for detoxification and 
metabohte excretion indeed increases, as expected by the up- 
regulation of transcripts associated with heme and ammonia 
detoxification/excretion (Figures 4-6), then the cytoskeletal 
dynamics may facilitate the movements of newly synthesized 
membrane-bound transporters that mediate the excretion of 
metabolites, such as ABC transporters. 

It is also possible that a more dynamic microtubule cytoskeleton 
would facilitate the movements of organelles within the epithelial 
cells. For example, retraction of mitocliroiidria from the apical 
microvilli in principal cells is associated with a decrease of fluid 
secretion in Malpighian tubules isolated from pupal stages of 
mosquitoes [85]. Thus, similar microtubule-mediated movements 
of mitochondria may occur during the chronic processing of blood 
meals (24-48 h after blood feeding) to further contribute to a 
putative decreased capacity for diuresis. 

Conclusions and significance 

The present study provides the first transcriptomic analysis of 
the Malpighian tubules of a mosquito after a blood meal, and is 
also the first to be conducted in A. albopictus after blood feeding. 
The results reveal molecular changes in transcript accumulation in 
the tubule epithelium within the first 24 h after a blood meal that 
suggest a remarkable functional transition of the epithelium from 
one dedicated to electrolyte and fluid excretion to one dedicated to 
detoxification and metabolite processing (Figure 8). Moreover, the 
results uncover new putative roles of the Malpighian tubules in the 
chronic processing of blood meals after the post-prandial diuresis 
ends ~2 h after a blood meal [8]. Thus, the tubule epithelium 
may represent an even more valuable target for the development 
of novel insecticides than has been previously appreciated. The 
next important step to complete is to validate the hypothesized 
functional transition of the epithelium after a blood meal using 
biochemical and physiological approaches. 

Supporting Information 

Figure SI Lists of transcripts that exhibited a sustained down- 
regulation after blood feeding and were associated with the 
'oxidative phosphorylation /ATP synthesis' DAVID functional 
cluster (see Table 2). Blue shading indicates significant down- 
regulation compared to NBF controls. Lack of shading indicates 
no significant difference relative to NBF controls. Degree of 
shading is based on value of Logg fold change as indicated by the 
scale below on the right. 
(TIF) 

Figure S2 Lists of transcripts that exhibited a sustained down- 
regulation after blood feeding and were associated with the 
'ATPase activity' DAVID functional cluster (see Table 2). Blue 
shading indicates significant down-regulation compared to NBF 
controls. Lack of shading indicates no significant difference 
relative to NBF controls. Degree of shading is based on value 
of Log2 fold change as indicated by the scale below on the 
right. 
(TIF) 



Figure S3 Lists of transcripts that exliibited a sustained down- 
regulation after blood feeding and were associated with the 
'glycolysis' DAVID functional cluster (see Table 2). Blue shading 
indicates significant down-regulation compared to NBF controls. 
Lack of shading indicates no significant difference relative to NBF 
controls. Degree of shading is based on value of Log2 fold change 
as indicated by the scale below on the right. 
(TIF) 

Figure S4 Lists of transcripts that exliibited a sustained down- 
regulation after blood feeding and were associated with the 'sugar/ 
inositol transporter' DAVID functional cluster (see Table 2). Blue 
shading indicates significant down-regulation compared to NBF 
controls. Lack of shading indicates no significant difference relative 
to NBF controls. Degree of shading is based on value of Logg fold 
change as indicated by the scale below on the right. 
(TIF) 

Figure S5 Lists of transcripts that exhibited a sustained down- 
regulation after blood feeding and were associated with the 
'protein biosynthesis' DAVID functional cluster (see Table 2). Blue 
shading indicates significant down-regulation compared to NBF 
controls. Lack of shading indicates no significant difference relative 
to NBF controls. Degree of shading is based on value of Logj fold 
change as indicated by the scale below on the right. 
(TIF) 

Figure S6 Lists of transcripts that exhibited a sustained down- 
regulation after blood feeding and were associated with the 
'translational elongation' DAVID functional cluster (see Table 2). 
Blue shading indicates significant down-regulation compared to 
NBF controls. Lack of shading indicates no significant difference 
relative to NBF controls. Degree of shading is based on value of 
Log9 fold change as indicated by the scale below on the right. 
(TIF) 

Figure S7 Lists of transcripts that exhibited a sustained up- 
regulation after blood feeding and were associated with the 
'thioredoxin' DAVID functional cluster (see Table 2). Yellow 
shading indicates significant up-regulation compared to NBF 
controls. Lack of shading indicates no significant difference relative 
to NBF controls. Degree of shading is based on value of Log2 fold 
change as indicated by the scale below on the right. 
(TIF) 

Figure S8 Lists of transcripts that exhibited a sustained up- 
regulation after blood feeding and were associated with the 
'glutathione S-transferase' DAVID functional cluster (see Table 2). 
Yellow shading indicates significant up-regulation compared to 
NBF controls. Lack of shading indicates no significant difference 
relative to NBF controls. Degree of shading is based on value of 
Log2 fold change as indicated by the scale below on the right. 
(TI^ 

Figure S9 Lists of transcripts that exhibited a sustained up- 
regulation after blood feeding and were associated with the 
'cofactor metabolic process' DAVID functional cluster (see 
Table 2). Yellow shading indicates significant up-regulation 
compared to NBF controls. Lack of shading indicates no 
significant difference relative to NBF controls. Degree of shading 
is based on value of Log2 fold change as indicated by the scale 
below on the right. 
(TIF) 

Figure SIO Lists of transcripts that exhibited a sustained up- 
regulation after blood feeding and were associated with the 
'vitamin biosynthetic process' DAVID functional cluster (see 
Table 2). Yellow shading indicates significant up-regulation 



PLOS Neglected Tropical Diseases | www.plosntds.org 



13 



June 2014 | Volume 8 | Issue 6 | e2929 



Changes to the Renal Transcriptome of Mosquitoes after a Blood Meal 



compared to NBF controls. Lack of shading indicates no 
significant difference relative to NBF controls. Degree of shading 
is based on value of Logg fold change as indicated by the scale 
below on the right. 
(TIF) 

Figure Sll Lists of transcripts that exhibited a sustained up- 
regulation after blood feeding and were associated with the 
'ATPase, AAA+ type' DAVID functional cluster (see Table 2). 
Yellow shading indicates significant up-regulation compared to 
NBF controls. Lack of shading indicates no significant difference 
relative to NBF controls. Degree of shading is based on value of 
Log2 fold change as indicated by the scale below on the right. 
(TIF) 

Figure SI 2 Lists of transcripts that exhibited a sustained up- 
regulation after blood feeding and were associated with the 
'proteasome complex' DAVID functional cluster (see Table 2). 
Yellow shading indicates significant up-regulation compared to 
NBF controls. Lack of shading indicates no significant difference 
relative to NBF controls. Degree of shading is based on value of 
Log2 fold change as indicated by the scale below on the right. 
(TIF) 

Figure S13 Lists of transcripts that exhibited a sustained up- 
regulation after blood feeding and were associated with the 'amine 
biosyntlietic process' DAVID functional cluster (see Table 2). 
Yellow shading indicates significant up-regulation compared to 
NBF controls. Lack of shading indicates no significant difference 
relative to NBF controls. Degree of shading is based on value of 
Log2 fold change as indicated by the scale below on the right. 
(TIF) 

Figure S14 Lists of transcripts that exhibited a sustained up- 
regulation after blood feeding and were associated with the 
'vitamin binding' DAVID functional cluster (see Table 2). Yellow 
shading indicates significant up-regulation compared to NBF 
controls. Lack of shading indicates no significant difference relative 

References 

1 . Benedict MQ, Lcvine RS, Hawley WA, Lounibos LP (2007) Spread of the tiger: 
global risk of inva.sion by the mosquito Aecfes albopictus. Veetor Borne Zoonotic 
Dis 7: 76-8,'j. 

2. Enserink M (2008) Entomology. A mosquito goes global. Science 320; 864— 
866. 

3. Roehlin I, Ninivaggi DV, Hutchinson ML, Farajollahi A (2013) Climate change 
and range expansion of the Asian tiger mosquito [Aedes albopictus) in Northeastern 
USA: implications for public health practitioners. PLoS One 8: e60874. 

4. Gratz NG (2004) Critical review of the vector status of Aedes albopictus. Med Vet 
Entomol 18: 21,5-227. 

5. VontasJ, Kioulos E, PavHdi N, Morou E, della Torre A, et al. (2012) Insecticide 
resistance in the major dengue vectors Aedes albopictus and Aedes aegrpti. Pesticide 
Biochemistry and Physiology 104: 126-131. 

6. Raphemot R, Rouhier MF, Hopkins CR, Gogliotu RD, Lovell KM, et al. (2013) 
Eliciting renal failure in mosquitoes with a small-molecule inhibitor of inward- 
rectilying potassium channels. PLoS One 8: e64905. 

7. Beyenbaeh KW, Piermarini PM (201 1) Transecllular and paraceilular pathways 
of transepithclial lluid secretion in Malpighian (renal) tubules of the yellow fever 
mosquito Aedes aegypti. Acta Physiol 202: 387—407. 

8. Williams JC, Hagedorn HH, Beyenbaeh KW (1983) Dynamic changes in flow 
rate and composition of urine during the post blood meal diuresis iji Aedes aegypti. 
J Comp Physiol [B] 1,53: 257-266. 

9. Isoe J, Scaraffia PY (2013) Urea Synthesis and Excretion in Aedes aegypti 
Mosquitoes Arc Regulated by a Unique Cross-Talk Mechanism. PLoS One 8: 
e65393. 

10. Briegel H, Lea AO (1975) Relationship between protein and proteolytic activity 
in the midgut of mosquitoes. J Insect Physiol 21; 1597—1604. 

11. Lemos FJ, Cornel AJ, Jacobs-Lorena M (1996) Trypsin and aminopeptidase 
gene expression is affected by age and food composition in Anopheles gambiae. 
Insect Biochem Mol Biol 26; 651-658. 

12. Bonizzoni M, Dunn WA, CampbeU CL, Olson KE, Dimon MT, et al. (201 1) 
RNA-seq analyses of blood-induced changes in gene expression in the mosquito 
vector species, Aedes aegypti. BMC Genomics 12: 82. 



to NBF controls. Degree of shading is based on value of Logj fold 

change as indicated by the scale below on the right. 

(TIF) 

Figure S15 Lists of transcripts that exhibited a sustained up- 
regulation after blood feeding and were associated with the 
'tubulin, GTPase domain' DAVID functional cluster (see Table 2). 
Yellow shading indicates significant up-regulation compared to 
NBF controls. Lack of shading indicates no significant difference 
relative to NBF controls. Degree of shading is based on value of 
Log2 fold change as indicated by the scale below on the right. 
(TIF) 

Table SI Number of 'raw reads' and 'reads mapped' for the 
Malpighian tubule cDNA libraries that were sequenced using 
RNA-Seq. 
(DOCX) 

Table S2 ANOVA of non-blood fed libraries across time points. 
Number of reads per transcript per library was used as the 
independent variable. 
(DOCX) 

Table S3 ANOVA of blood fed libraries across time points. 
Number of reads per transcript per library was used as the 
independent variable. 
(DOCX) 

Acknowledgments 

Wc thank Dr. Asela Wijcratne, Dr. Saraiiga Wijeratiie, Dr. Tea Meulia, 
Dr. Matthew F. Rouhier, and Ms. Nuris Acosta (all from the Ohio State 
University) for their assistance and guidance with this project. 

Author Contributions 

Conceived and designed the experiments: PMP CJE BJC. Performed the 
experiments: CJE. Analyzed the data: PMP CJE BJC. Contributed 
reagents/materiais/anafysis tools: PMP CJE BJC. Wrote tire paper: PMP 
CJE BJC. 



13. Marinotti O, Calvo E, Nguyen QK, Dissanayakc S, Ribciro JM, ct al. (2006) 
Gcnomc-wide analysis of gene expression in adult Anopheles gambiae. Insect Mol 
Biol 15: 1-12. 

14. Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, et al. (2002) 
The genome sequence of tlie malaria mosquito Anopheles gambiae. Science 298: 
129-149. 

15. Ribeii'oJM (2003) A catalogue oi' Anopheles gambiae transcripts signilicantly more 
or less expressed following a blood meal. Insect Biochem Mol Biol 33: 865-882- 

16. Dana AN, Hong YS, Kern MK, HiUenmcycr ME, Harkcr BW, et al. (2005) 
Gene expression patterns associated with blood-feeding in the malaria mosquito 
Anopheles gambiae. BMC Genomics 6: 5. 

1 7. Price DP, Nagarajan V, Churbanov A, Houde P, Milligan B, et al. (201 1) The 
fat body transcriptomcs of the yellow fever mosquito Aedes aegypti, pre- and post- 
blood meal. PLoS One 6: e22573. 

18. Rinkcr DC, Pitts RJ, Zhou X, Suh E, Rokas A, et al. (2013) Blood meal-induced 
changes to antcnnal transcriptome profiles reveal shifts in odor sensitivities in 
Anopheles gambiae. Proc Nad Acad Sci U S A 110: 8260-8265. 

19. Sanders HR, Evans AM, Ross LS, Gill SS (2003) Blood meal induces global 
changes in midgut gene expression in the disease vector, Aedes aegppti. Insect 
Biochem Mol Biol 33: 1105-1122. 

20. Thangamani S, Wikel SK (2009) Differential expression ot~ Aedes aegypti salivary 
transcriptome upon blood feeding. Parasit Vectors 2: 34. 

21. Pannabecker TL, Hayes TK, Beyenbaeh KW (1993) Regulation of epitliehal 
shunt conductance by the peptide leucokinin. J Membr Biol 132: 63-76. 

22. Drake LL, Boudko DY, Marinotti O, Carpenter VK, Dawe AL, ct al. (2010) 
The Aquaporin gene family of the yellow fever mosquito, Aedes aegypti. PloS one 
5: cl5578. 

23. Chomczynski P, Sacchi N (1987) Single-step method of RNA isolation by acid 
guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem 162: 
156-159. 

24. GoecksJ, NekiTitcnko A, Taylor J (2010) Galaxy: a comprehensive approach for 
supporting accessible, reproducible, and transparent computational research in 
the life sciences. Genome Biol 11: R86. 



PLOS Neglected Tropical Diseases | www.plosntds.org 



14 



June 2014 | Volume 8 | Issue 6 | e2929 



Changes to the Renal Transcriptome of Mosquitoes after a Blood Meal 



25. Martin M (2011) Cutadapt removes adapter sequences from high-throughput 
sequencing reads. EMBnetJ 17: 10-12. 

26. Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows- 
Wheeler transform. Bioinformatics 25: 1754—1760. 

27. Nene V, Wortman JR, Lawson D, Haas B, Kodira C, et al. (2007) Genome 
sequence of Aedes aegypti, a major arbovirus vector. Science 316: 1718—1723. 

28. Poelchau MF, Reynolds JA, Denlinger DL, Elsik CG, Armbruster PA (201 1) A 
de novo transcriptome of tlic Asian tiger mosquito, Aedes albopictus, to identify 
candidate transcripts for diapause preparation. BMC Genomics 12: 619. 

29. Bonizzoni M, Dunn WA, CampbeU CL, Olson KE, Marinotti O, et al. (2012) 
Strain Variation in the Transcriptome of the Dengue Fever Vector, Aedes aeg)pti. 
G3 (Bethesda) 2: 103-114. 

30. Akbari OS, Antoshechkin I, Amrhein H, Williams B, Diloreto R, et al. (2013) 
The developmental transcriptome of the mosquito Aedes aegypti, an invasive 
species and major arbovirus vector. G3 (Bethesda) 3: 1493—1509. 

31. Anders S, Hubcr W (2010) Differential expression analysis for sequence count 
data. Genome Biol 11: R106. 

32. Huang da W, Sherman BT, Lcmpicki RA (2009) Systematic and integrative 
analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4: 
44-57. 

33. Huang da W, Sherman BT, Lcmpicki RA (2009) Bioinformatics enrichment 
tools: paths toward the comprehensive functional analysis of large gene lists. 
Nucleic Acids Res 37: 1-13. 

34. Beyenbach KW, Pannabecker TL, Nagel W (2000) Central role of the apical 
membrane H^-ATPase in electrogenesis and epithelial transport in Malpighian 
tubules. J Exp Biol 203: 1459-1468. 

35. Weng XH, Huss M, Wieczorek H, Beyenbach KW (2003) The V-type H""- 
ATPase in Malpighian tubules of Aedes aegypti: localization and activity. J Exp 
Biol 206: 2211-2219. 

36. Beyenbach KW (2003) Transport mechanisms of diuresis in Malpighian tubules 
of insects. J Exp Biol 206: 3845-3856. 

37. Beyenbach KW, Wieczorek H (2006) The V-type H^ ATPase: molecular 
structure and function, physiological roles and regulation. J Exp Biol 209: 577- 
589. 

38. Wu DS, Beyenbach KW (2003) The dependence of electrical transport pathways 
in Malpighian tubules on ATP. J Exp Biol 206: 233-243. 

39. Scott BN, Yu MJ, Lee LW, Beyenbach KW (2004) Mechanisms of K+ transport 
across basolateral membranes of principal cells in Malpighian tubules of the 
yellow fever mosquito, Aedes aegypti. J Exp Biol 207: 1655-1663. 

40. Beyenbach KW, Masia R (2002) Membrane conductances of principal cells in 
Malpighian tubules oi' Aedes aegypti.J Insect Physiol 48: 375-386. 

41. Piermarini PM, Rouhier MF, Schepel M, Kosse C, Beyenbach KW (2013) 
Cloning and functional characterization of inward-rectifying potassium (Kir) 
channels fi'om Malpighian tubules of the mosquito Aedes aegypti. Insect Biochem 
Mol Biol 43: 75-90. 

42. Rodan AR, Baum M, Huang CL (2012) The Drosophila NKCC Ncc69 is 
required for normal renal tubule function. Am J Physiol Cell Physiol 303: C883- 
894. 

43. WilUamsJC, Beyenbach KW (1984) Dilferential clfects of secretagogues on the 
electrophysiology of the Malpighian tubules of tlie yellow fever mosquito. J Comp 
Physiol [B] 154: 301-309. 

44. Sawyer DB, Beyenbach KW (1985) Dibutyryl-cAMP increases basolateral 
sodium conductance of mosquito Malpighian tubules. AmJ Physiol Regul Integr 
Comp Physiol 248: R339-R345. 

45. PuUikuth AK, Aimanova K, Kang'ethe W, Sanders HR, Gill SS (2006) 
Molecular characterization of sodium/proton exchanger 3 (NHE3) from the 
yellow fever vector, Aedes aeg)ipti.] Exp Biol 209: 3529-3544. 

46. Liu L, Johnson WA, Welsh MJ (2003) Drosophila DEG/ENaC pickpocket genes 
are expressed in the tracheal system, where they may he involved in liquid 
clearance. Proc Nad Acad Sci U S A 100: 2128-2133. 

47. Cameron P, Hiroi M, NgaiJ, Scott K (2010) The molecular basis for water taste 
in Drosophila. Nature 465: 91-95. 

48. Piermarini PM, Grogan LF, Lau K, Wang L, Beyenbach KW (2010) A SLC4- 
like anion exchanger from renal tubules of the mosquito [Aedes aegypti): evidence 
for a novel role of stellate cells in diuretic fluid secretion. Am J Physiol Regul 
Integi- Comp Physiol 298: R642-660. 

49. Linser PJ, Neira Oviedo M, Hirata T, Seron TJ, Smith KE, et al. (2012) Slc4- 
like anion transporters of the larval mosquito alimentary canal. J Insect Physiol 
58: 551-562. 

50. Liu K, Tsujimoto H, Cha SJ, Agre P, Rasgon JL (2011) Aquaporin water 
channel AgAQPl in the malaria vector mosquito Anopheles gambiae during blood 
feeding and humidity adaptation. Proc Natl Acad Sci U S A 108: 6062-6066. 

51. Tsujimoto H, Liu K, Linser PJ, Agre P, Rasgon JL (2013) Organ-specific splice 
variants of aquaporin water channel AgAQPl in the malaria vector Anopheles 
gambiae. PLoS One 8: e75888. 

52. Niu LL, Fallon AM (2000) Differential regulation of ribosomal protein gene 
expression in Aedes aegypti mosquitoes before and after the blood meal. Insect Mol 
Biol 9: 613-623. 

53. Graca-Souza AV, Maya-Monteiro C, Paiva-Silva GO, Braz GR, Paes MC, et al. 
(2006) Adaptations against heme toxicity in blood-feeding arthropods. Insect 
Biochem Mol Biol 36: 322-335. 

54. Pascoa V, Oliveira PL, Dansa-Petretski M, SilvaJR, Alvarenga PH, et al. (2002) 
Aedes aegypti peritrophic matrix and its interaction with heme during blood 
digestion. Insect Biochem Mol Biol 32: 517-523. 



55. Devenport M, Alvarenga PH, Shao L, Fujioka H, Bianconi ML, et al. (2006) 
Identification of the Aedes aegfpli Peritrophic Matrix Protein AelMUCI as a 
Heme-Binding Protein. Biochemistry 45: 9540-9549. 

56. Pcreii-a LO, Oliveira PL, Almeida IC, Paiva-Silva GO (2007) Biglutaminyl- 
biliverdin IX alpha as a heme degradation product in the dengue fever insect- 
vector Aedes aegypti. Biochemistry 46: 6822-6829. 

57. Han Beerntsen BT, Li J (2007) The tryptophan oxidation pathway in 
mosquitoes with emphasis on xanthurenic acid biosynthesis. J Insect Physiol 53: 
254-263. 

58. Lima VL, Dias F, Nunes RD, Percira LO, Santos TS, et al. (2012) The 
antioxidant role of xanthurenic acid in the Aedes aegypti midgut during digestion 
of a blood meal. PLoS One 7: e38349. 

59. Lumjuan N, Stevenson BJ, Prapanthadara LA, Somboon P, Brophy PM, et al. 
(2007) The Aedes aegypti glutathione transferase family. Insect Biochem Mol Biol 
37: 1026-1035. 

60. Toh SQ, Glanfield A, Gobert GN, Jones MK (2010) Heme and blood-feeding 
parasites: friends or foes? Parasit Vectors 3: 108. 

61. Dermauw W, Van Leeuwen T (2014) The ABC gene family in arthropods: 
Comparative genomics and role in insecticide transport and resistance. Insect 
Biochemistry and Molecular Biology 45: 89—1 10. 

62. Bariami V, Jones CM, Poupardin R, Vontas J, Ranson H (2012) Gene 
amplification, ABC transporters and cytochrome P450s: unraveling the 
molecular basis of pyrethroid resistance in the dengue vector, Aedes aegipti. 
PLoS Negl Trop Dis 6: el692. 

63. Strode C, Wondji CS, David JP, Hawkes NJ, Lumjuan N, et al. (2008) Genomic 
analysis of detoxification genes in the mosquito Aedes aegypti. Insect Biochem Mol 
Biol 38: 113-123. 

64. Boyer TD, Vessey DA, Holcomb C, Saley N (1984) Studies of the relationship 
l^etween die catalytic activity and binding of non-sul:)stratc ligands by the 
glutathione S-transferases. BiochcmJ 217: 179-185. 

65. Phillips O, Mande TJ, Tuffery AR, Hcyworth CM, Wilson SR, et al. (1984) On 
the possible role of biliverdin stimulation of cycHc AMP levels as a trigger for 
liver regeneration in the rat. Biochemical Pharmacology 33: 1963-1967. 

66. Caccuri AM, Aceto A, Piemonte F, Di Ilio C, Rosato N, et al. (1990) Interaction 
of hemin with placental glutathione transferase. Eur J Biochem 189: 493-497. 

67. Cole SPC, Deeley RG (2006) Transport of glutathione and glutathione 
conjugates by MRPl. Trends in Pharmacological Sciences 27: 438—446. 

68. Krishnamurthy P, Xie T, Schuetz JD (2007) The role of transporters in ceUular 
heme and porphyrin homeostasis. Pharmacology & Therapeutics 1 14: 345—358. 

69. Basseville A, Bates SE (201 1) Gout, genetics and ABC transporters. FIOOO Biol 
Rep 3: 23. 

70. Wang L, Kiuchi T, Fujii T, Daimon T, Li M, et al. (2013) Mutation of a novel 
ABC transporter gene is responsible for the failure to incorporate uric acid in the 
epidermis of ok mutants of die silkworm, Bonibyx mori. Insect Biochemistry and 
Molecular Biology 43: 562-571. 

71. Tarling EJ, VaUim TQdA, Edwards PA (2013) Role of ABC transporters in lipid 
transport and human disease. Trends in endocrinology and metabolism: TEM 
24: 342-350. 

72. Bellarosa C, Bortolussi G, Tiribelh C (2009) The Role of ABC Transporters in 
Protecting Cells from Bilirubin Toxicity. Current Pharmaceutical Design 15: 
2884-2892. 

73. Suraweera A, Munch G, Hanssum A, Bcrtolotti A (2012) Failure of amino acid 
homeostasis causes cell death following proteasome inhibition. Mol Cell 48: 242— 
253. 

74. Scaraffia PY, Isoe J, MurUlo A, Wells MA (2005) Ammonia metabolism in Aedes 
aegypti. Insect Biochem Mol Biol 35: 491-503. 

75. Scaraffia PY, Zhang Wysocki VH, Isoe J, Wells MA (2006) Analysis of whole 
body ammonia metabolism in Aedes aegypti using [^"^N]-labeled compounds and 
mass spectrometry. Insect Biochem Mol Biol 36: 614—622. 

76. Scaraffia PY, Tan G, Isoe J, Wysocki VH, Wells MA, et al. (2008) Discovery of 
an alternate metabolic pathway for urea synthesis in adult Aedes aegypti 
mosquitoes. Proc Nati Acad Sci U S A 105: 518-523. 

77. Scaraffia PY, Zhang Thorson K, Wysocki VH, Miesfeld RE (2010) 
Differential ammonia metabolism in Aedes aegypti fat body and midgut tissues. 
J Insect Physiol 56: 1040-1049. 

78. Scaraffia PY, WcUs MA (2003) Proline can be utilized as an energy substrate 
during flight of Aedes aegypti females. J Insect Physiol 49: 591—601. 

79. Briegcl H (1986) Protein catabolism and nitrogen partitioning during oogenesis 
in the mosquito Aedes aegypti. Journal of Insect Physiology 32: 455—462. 

80. Bursell E (1965) Nitrogenous waste products of the tsetse fly, Glossina morsitans. 
J Insect Physiol 11: 993-1001. 

81. Wigglesworth VB (1931) The Physiology Of Excretion In A Blood-Sucking 
Insect, Rhodnius Prolixus (Hemiptera, Reduviidae): I. Composition Of The Urine. 
Journal of Experimental Biology 8: 41 1-427. 

82. O'Donnell MJ, MaddrcU SHP, Gardiner BOC (1983) Transport of uric acid by 
the Malpighian tubules of Rhodnius prolixus and other insects. Journal of 
Experimental Biology 103: 169-184. 

83. Horgan CP, McCaffrey MW (2011) Rab GTPases and microtubule motors. 
Biochem Soc Trans 39: 1202-1206. 

84. Karas K, Brauer P, Petzel D (2005) Actin redistribution in mosquito Malpighian 
tubules after a blood meal and cyclic AMP stimulation. J Insect Physiol 51: 
1041-1054. 

85. Bradley TJ, Snyder C (1989) Fluid secretion and microvillar ultrastructure in 
mosquito Malpighian tubules. AmJ Physiol 257: R1096-1 102. 



PLOS Neglected Tropical Diseases | www.plosntds.org 



15 



June 2014 | Volume 8 | Issue 6 | e2929 



Changes to the Renal Transcriptome of Mosquitoes after a Blood Meal 



86. O'Connor KR, Beyenbach KW (2001) Chloride channels in apical membrane 
patches of stellate cells of Malpighian tubules of Aetfes aeg}>pti. J Exp Biol 204; 
367-378. 

87. Piermarini PM, Hinc RM, Schcpcl M, Miyauchi JT, Beyenbach KW (2011) 
Role of an apical K,C1 cotransporter in urine formation by renal tubules of the 
yellow fever mosquito {Asdes aegypti). Am J Physiol Ret^l Integr Comp Physiol 
301: R1318-1337. 



88. Xiang MA, Linser PJ, Price DA, Harvey WR (20 1 2) Locahzation of two Na+- or 
K^-H^ antiporters, AgNHAl and AgNHA2, in Anopheles gamhme- larval 
Malpighian tubules and the functional expression of AgNHA2 in yeast. J Insect 
Physiol 58: 570-579. 

89. HegartyJL, Zhang B, Carroll MC, Cragoe EJJ, Beyenbach KW (1992) Effects of 
amiloridc on isolated Malpighian tubules of the yellow fever mosquito {Aedes 
aeg)pi].J Insect Physiol 38: 329-337. 



PLOS Neglected Tropical Diseases | www.plosntds.org 



16 



June 2014 | Volume 8 | Issue 6 | e2929 



